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Heated high speed subsonic and supersonic jets operating on- or off-design are a source 
of noise that is not yet fully understood. Helium-air mixtures can be used in the correct 
ratio to simulate the total temperature ratio of heated air jets and hence have the potential 
to provide inexpensive and reliable flow and acoustic measurements. This study presents a 
combination of flow measurements of helium-air high speed jets and numerical simulations 
of similar helium-air mixture and heated air jets. Jets issuing from axisymmetric conver- 
gent and convergent-divergent nozzles are investigated, and the results show very strong 
similarity with heated air jet measurements found in the literature. This demonstrates the 
validity of simulating heated high speed jets with helium-air in the laboratory, together 
with the excellent agreement obtained in the presented data between the numerical predic- 
tions and the experiments. The very close match between the numerical and experimental 
data also validates the frozen chemistry model used in the numerical simulation. 


Nomenclature 


a Speed of sound 

c Specific heat 

D Nozzle exit diameter 

H Enthalpy 

M Mach number 

M a Acoustic Mach number ( Uj / a^) 

Md Nozzle design Mach number 

P Pressure 

P p Pitot pressure 

R Gas constant 

r Radial direction, coordinate 

S Entropy 

T Temperature 

TTR Total temperature ratio (T t / TV) 

U Mean velocity 

x Downstream direction, coordinate 

Greek 

Mass fraction 
p Density 

X Molar fraction 

Subscript 
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t Post-doctoral researcher, Aerospace Engineering Department. 
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1 Local quantity 

2 Quantity measured behind the shock created by the pitot probe 
j Jet exit plane mean quantity 

t Total or stagnation quantity 
oo Atmospheric quantity 

Superscript 

air Quantity relating to air 
He Quantity relating to helium 
heated Quantity relating to heated air 
mix Quantity relating to a helium-air mixture 


I. Introduction 

Military aircraft jet engines have noise characteristics much louder than civilian aircraft due to their very 
low bypass ratio, high exhaust temperature, and high exhaust velocities. The high intensity noise generated 
by these heated subsonic or supersonic jets is a health hazard to ground crews as well as an annoyance to 
communities in the vicinity of airfields. This has led to a need for reduction of jet noise by developing noise 
suppression mechanisms that include new nozzle design concepts such as chevrons, corrugations, beveled 
nozzles, or other non-axisymmetric geometries. Research by Tanna et a/., 1 Lau, 2 Seiner et a/., 3 and Tam 
and Chen 4 demonstrated that there are significant differences in how moderately heated supersonic jets 
generate noise compared with their unheated counterparts. Mixing noise is present in all jets and is caused 
by the coherent and incoherent turbulent structures. Adding some heating to a supersonic air jet helps 
to stabilize the instability waves in the shear layer, but when these instability waves travel faster than 
the ambient speed of sound, Mach wave radiation occurs in addition to the mixing noise. Jets operating 
supersonically and off-design also exhibit shock noise which is dominant in the upstream direction. 

To study jets under realistic engine operating conditions, they should be heated so that the convection 
velocity of the turbulent structures is supersonic and the effects of Mach wave radiation can be measured. 
In practical terms, a heated jet facility involves significant complexity and expenses. To avoid this, the low 
density and high velocity of hot jets have been simulated in previous studies 5-7 by using a lower density gas 
with different properties than air, namely, helium. By using a helium-air mixture, it becomes possible to 
simulate either the density or the acoustic velocity of a heated jet without having to heat the air. Although 
there are some differences between an actual heated jet and the above described simulated experiment, the 
major features of the noise generation process seem well represented and provided very satisfactory results 
in previous studies. 5-7 However, as noise reduction concepts are being investigated, with increasing demand 
on the accuracy of the acoustic and flow measurements, a detailed assessment of the limits of the simulation 
of hot air jets with helium-air mixtures needs to be made. The present study investigates the details of the 
flow properties of both heated and heat simulated jets, using a combination of experimental measurements 
and Computational Fluid Dynamic (CFD) simulations. Experiments performed with cold, pure air jets and 
with helium- air mixtures highlight the effects of simulated heat on the jet flow properties. Numerical results 
with cold and heated pure air jets as well as helium-air mixture jets are used to assess the validity of the 
helium addition to properly simulate heat as well as validate a multi-species diffusion model. 

II. Approach 

The purpose of this study is to provide detailed comparisons between the flow properties of helium-air 
mixture jets simulating heated jets with heated jets. A short description of the simulation of heated jets 
with a helium-air mixture is presented. Then the experimental and computational approaches are described. 

A. Simulating Heating with Helium- Air Mixtures 

In order to make acoustic measurements that can be directly compared to aircraft engine measurements, 
the temperature of the jet is an important parameter that needs to be replicated. A heated jet has different 
physical characteristics than a cold jet, due to the increase in jet exit velocity and decrease in jet density. 
This affects the acoustical properties of the jet. In order to avoid the power, infrastructure, and operating 
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cost requirements of a heated jet facility, a helium-air mixture can be used as a means to simulate the desired 
jet properties without actually heating the air flow. 8 This methodology was developed at the Pennsylvania 
State University by Kinzie and McLaughlin 6 and Doty and McLaughlin. 5 A similar methodology is being 
used by Papamoschou 7 at the University of California Irvine. The quantities that need to be matched when 
doing acoustic measurements in hot jets are the jet acoustic velocity aj and the jet density pj. Helium was 
chosen for its low density and high acoustic velocity and the fact that it is an inert gas which is inexpensive 
to obtain and handle. In order to account for varying atmospheric conditions, the characterizing parameter 
is the jet exit to ambient temperature ratio Tj/T a . For the helium-air jet to match a corresponding heated 
jet, two non-dimensional groups are required: the density ratio pj/p a and the acoustic velocity ratio aj/a a . 
It is possible to match both, as mentioned by Kinzie, 6 by mixing some helium with ambient air. However, it 
is impractical in a large facility. In prior experimental work by Doty and McLaughlin, 5 the two parameter 
matching methods were investigated, with direct comparisons to heated air measurements. Only very small 
variations in the acoustic spectra were observed and were of the same order of magnitude as the experimental 
uncertainty of the measurements. More recent careful comparisons 9 with measurements performed in other 
facilities have shown very good agreement when matching the acoustic velocity of the mixture jet to that of 
a heated jet. Therefore, matching of the acoustic velocity is chosen for the experimental and numerical work 
presented in this study, meaning the mixture properties are defined such that: 

/ \ heated / \ mix 

= ° V~f RT 3 = \J ^ mix R mix Tj nix (1) 

More details on the effectiveness of helium in simulating hot jets, as well as the implications associated 
with this methodology are discussed in Doty and McLaughlin 5 and Papamoschou. 7 These two references 
discuss the safety and economic benefits of helium-air mixtures compared to heating air for the experiments 
conducted in this facility. 

B. Experimental Approach 

Experimental measurements of the mean flow properties were performed using Pitot probes within the jet 
plume. The schematic diagram presented in Fig. la shows the nomenclature used for the flow quantities 
around a single probe. A subscript ( )2 indicates the properties just downstream of the normal shock 
portion of the shock wave, immediately upstream of the Pitot probe, and the subscript ( )i indicates the 
local properties ahead of the shock. Figure lb and lc show photographs of the Pitot probes used for the 
measurements. It consists of a rake of five probes of outer diameter 0.6mm (0.025”) and inner diameter 
0.25mm (0.01”). Apart from this new probe rake, the overall experimental setup is the same as the one 
described in Miller et al. 10 


*0 p 

M 

Pa 

M x 






Figure 1. a) Schematic diagram of a Pitot pressure probe in a free stream flow, b) detailed view of the Pitot 
rake at the exit plane of a round nozzle, c) Pitot setup. 
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The Pitot pressure P p , can be used to find a value for the local Mach number M\ using a similar process 
developed for pure air jets. Wherever the flow is supersonic, a shock forms in front of the probe, and P p is 
equal to the total pressure behind the shock, P t 2. Mi is then calculated using the Rayleigh Pitot formula 
shown in Eq. (2) written with static pressure before the shock. Wherever the flow is subsonic, there is no 
shock ahead of the probe, thus P p is equal to P t i and the isentropic flow formula of Eq. (3) can be used. 
The jets evaluated in this study are fully-expanded, the static pressure Pi can be assumed to be constant 
and equal to the ambient static pressure P ^ everywhere within the jet. A more detailed description of 
this processing is presented in Miller et al 10 as well as some detailed comparisons between simulations and 
experimental measurements. 
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However, attempting to perform this calculation in a helium-air mixture jet leads to a new complication: 
the mass concentration of helium ^ e is not constant within the mixing layer, leading to fluctuations of 
properties such as the local specific heat ratio 71 across the jet. Since no measurement of the species 
concentration is available within such high speed jets, an analytical approach needs to be used in order to 
yield estimates of the gas constitution. In steady axisymmetric turbulent jets the conservation of momentum 
in the x (streamwise) direction can be expressed as follow: 


du du d f du\ dP 
pU frc +pV lfr - dr ) = 


(4) 


where p t represents the eddy viscosity, or the effective viscosity of a turbulent fluid, and u and v are 
respectively the axial and radial components of the velocity. Note that the normal stress contribution is 
neglected. In the free-stream jets considered, the streamwise gradient of pressure can be neglected (dP/dx), 
such that the right hand side term is dropped. Then, within the jet mixing layer, three different phenomena 
can cause mixing of the jet inner gas mixture with the ambient air: turbulent mixing, molecular diffusion, 
and thermal diffusivity. Turbulent mixing can be described, as shown by Spalding, 11 using the conservation 
of mass of the entrained flow. The resulting equation relates the mass concentration of helium & He to the 
spatial velocity field: 


pu- 


d$ He 

dx 


+ pv- 


g$ He 

dr 


d_ 

dr 


p t d& He \ 
Set dr ) 


(5) 


where Set refers to the turbulence Schmidt number. When mixing occurs via molecular and thermal diffusion 
(of helium), rather than turbulent mixing, then the following equation, based on Fick’s law and derived by 
Schlichting 12 is used: 
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where D 12 denotes the coefficient of binary diffusion of helium in air and kr the thermal diffusion ratio. 
However, even though the jets of interest are not isothermal, they are unheated and therefore do not have 
large enough temperature gradients to produce appreciable thermal diffusion relative to other mechanisms. 
The thermal diffusion term can therefore be neglected. The remaining molecular diffusivity term is compared 
with the right hand side term of Eq. (5). The coefficient of molecular diffusion of helium in air is a known 
quantity, with value D 12— 62.10 -6 m 2 / s at ambient conditions. Therefore, the term pD 12 is of the order of 
magnitude 10 -4 . Also, the turbulence kinematic viscosity v t can be estimated, as shown by Schlichting 12 
and by Spalding 11 as iy t =0.026xro^xU c: where ro.5 is the half velocity point and U c is the centerline velocity 
of the jet. This leads to a dynamic viscosity p t with an order of magnitude 10 -1 for the jets considered 
(. D = 0.0127 m and Uj = 300 to 600 m/s). Using an estimated value of 0.7 for the turbulent Schmidt 
number, Sc tj as suggested by Panchapakesan and Lumley, 13 the constant in the right hand side of Eq. (5) 
has an order of magnitude of 10“ 1 . This is much higher than the order of magnitude 10 -4 obtained for the 
molecular diffusion. Therefore, molecular diffusion is neglected, which leaves turbulent mixing as the sole 
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contributor to the diffusion of helium outside of the jet. Eq. (4) and Eq. (5) are therefore the only equations 
to consider. Previous studies 13,14 have shown that while the value of Sc t varies across the shear layer, a 
constant value of 0.7 can be used and produces adequate concentration estimates. Therefore, Sc t and /i t are 
assumed constant and moved outside the partial derivatives. Furthermore, assuming the mass concentration 
of helium can be expressed as a sole function of the local flow velocity Lb, the chain rule for partial derivative 
can be applied to Eq. (5), and then combined with Eq. (4). These mathematical operations yield: 

( 1 \d<f> He d 2 u ! d 2 ^ He fdu i y == 

\Sc t J dUi dr 2 + OUt V dr ) U 

In order to obtain a solution for this second order differential equation the function U\(r) must be known. 
It was observed in past experimental studies 15 that the velocity profile in a jet follows a Gortler error function 
distribution, with amplitude proportional to the jet velocity Uj. The first and second derivative of such a 
function is also proportional to Uj. As a result, the second term of Eq. (7) is proportional to the square 
of the convection velocity, which means the first term can be neglected. Eq. (7) then reduces to a simple 
linear variation of the mass concentration of helium with respect to velocity. Following this development, 
reduction of the pressure data is performed with the underlying assumption that the helium concentration 
varies linearly with local velocity. Numerical and experimental results shown in this paper will illustrate the 
validity of this argument. 

The processing of the experimental data is performed with Matlab. The pitot pressure P p is first used for 
(using the Rayleigh and isentropic formula in the appropriate parts of the jet) an estimated value of the local 
Mach number M\ by using 7 j and Tj everywhere within the jet. Then, a first estimate of the local velocity 
is obtained, as well as the local static temperature, T\ using the Crocco relation. A better estimate of the 
velocity is then computed using T\. Then the mass and molar concentration of the species are evaluated, 
leading to an estimate for the local specific heat ratio and gas constant of the mixture. A refined velocity is 
computed from these values and the whole process is repeated until it converges to the final values for the 
local velocity U\ and temperature T\. Any other local property of the flow can be computed from these. 
The overall processing methodology is summarized in the flow chart of Fig. 2. 

C. Computational Approach 

In the previous section the experimental method and approach to find field- variables from helum-air jets is 
described. A numerical approach to find the same field-variables at the same spatial locations is discussed 
here. This numerical approach uses the Reynolds Averaged Navier-Stokes (RANS) equations, which are 
numerically integrated until a steady solution is achieved with the NPARC Alliance Wind-US 2.0 solver. 
Wind-US was developed by the NPARC Alliance which is a parternship between NASA Glenn Research 
Center and the U. S. Air Force Arnold Engineering Development Center, and additional contributors. A 
large number of validation cases of a broad range of flow-regimes have been made with the Wind-US solver 
by the NPARC Alliance. Wind-US can solve the Euler or Navier-Stokes equations in conjunction with many 
different turbulence models on structured or unstructured multi-block domains. It also has the ability to 
perform the computations in serial or parallel using the Message Passing Interface or the Parallel Virtual 
Machine. The Menter 16 Shear Stress Transport (SST) turbulence model is used in the present simulations 
for both heated air and helium-air mixtures to close the steady RANS equations. In addition, transport 
equations are solved for each species of the flow. In this case a helium, oxygen, and nitrogen mixture is 
chosen to represent the helium-air mixture. The traditional ideal gas model is also used for corresponding 
heated air simulations with a constant 7. The turbulent Prandtl number is set in all simulations at 0.90 and 
Sc t is set as 0.70. 

Generally CFD simulations of nozzles use ideal gas such as air because corresponding laboratory mea- 
surements use unheated or heated air. CFD simulations assume that the gas constant and ratio of specific 
heats are constant. Unfortunately, this assumption does not allow for the variation of the thermodynamic 
quantities inside and outside the nozzle, where they differ because of the large temperature differences. In 
aircraft jet engines the combustion process adds additional species to air changing the thermodynamic quan- 
tities. Simulation of these combustion reactions is computationally very expensive. However, in conventional 
aircraft engines the majority of the combustion process is complete when the gas reaches the nozzle inlet. 
Therefore, a frozen chemistry model can be used. The assumption can also be made, as discussed in the 
previous section, that diffusion is the only mechanism that changes the species concentration. Similar as- 


5 of 22 


American Institute of Aeronautics and Astronautics 



Rayleigh formula 



Figure 2. Flowchart of the pitot measurement processing. 


sumptions were made by Yoder et al . 17 with comparison between the effects of using a ideal gas (air), ideal 
gas mixtures, and a frozen exhaust gas mixture of subsonic and supersonic jets. 

In a helium-air mixture jet there is no combustion process and the number of species are relatively small 
compared to those produced in an aircraft engine. Since the helium-air mixture is known in the experiments 
and the only mixing occurs upstream of the nozzle, the mass species concentration is well known at the 
nozzle outlet. Far outside of the jet, the gas mixture corresponds to ambient air. Since no chemical process 
occurs in the shear layer of the jet, and mixing is mainly produced by turbulence, use of frozen chemistry in 
Wind-US is used. Frozen chemistry essentially assumes that the reaction source terms are zero and in these 
simulations only the diffusion and turbulent mixing of individual species is modeled. 

Wind-US has the ability to perform reacting flow calculations with a preexisting chemistry model and is 
initialized with the addition of a keyword in the Wind-US input file. Frozen chemistry modeling was used 
for the helium-air simulations presented in this study, but not for the heated air simulations. Additional 
inputs required by the chemistry module consist of the species concentration in the far-held and at the 
inflow boundary at the nozzle inlet. The full implementation of the frozen chemistry module and associated 
transport equations without reaction sources can be found in the Wind-US Documentation 18 or the Wind-US 
Developer’s Reference. 

Wind-US has a number of standard chemistry packages that contain various models of standard molecules. 
Unfortunately, none of these packages contain information for the element helium. An external chemistry hie 
is developed by Joo Kyung Suh of the Raytheon Corporation that contains the species helium along with the 
two most abundant elements in air, oxygen and nitrogen. The chemistry hie contains the thermodynamic 
coefficients of the three species, derived from McBride et al . 19 where three curves are dehned for each species. 
Each of these species has a set of three curves that describe the specihc heat at constant pressure, enthalpy, 
and entropy over a given range of temperatures in terms of coefficients cq and b{. Ranges of temperature for 
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the three curves are from 200 K to 20,000 K. Thermodynamic values defined by McBride et al. 19 based on 
the coefficients in the helium-air chemistry file are, 

Cp/ R = cl\T 2 + a 2 T 1 + <23 + U4T + a^T 2 + clqT^ + u 7 T 4 + <2 8 T 5 (8) 

tf 0 / (iET) = —a\T~ 2 + a 2 T _1 InT + a 3 + a 4 T/2 + a 5 T 2 /3 + a 6 T 3 /4 + a 7 T 4 /5 + a 8 T 5 /6 + b 1 /T (9) 

S%R = b 2 - aiT~ 2 /2 - a 2 T _1 + a 3 InT + a 4 T + a 5 T 2 /2 + a 6 T 3 / 3 + a 7 T 4 /4 + a 8 T 5 /5 (10) 

where c p is the specific heat at constant pressure, H is the enthalpy, R is the universal gas constant, S is 
the entropy, and T is the temperature. The chemistry file also contains the reaction rate coefficients for 
the three species with which Wind-US calculates the forward and backward reaction rates. Both rates are 
independent and are of the form, 


k = CT S exp [~D/K b T] (11) 

where K 5 is the Boltzmann constant, k is either the forward or backward reaction rate, and C and D are 
reaction rate coefficients. Since only a three species model is available, other species that are typically present 
in small quantity in air, such as Argon, are ignored. A similar assumption was made successfully by Yoder et 
al. 17 The frozen chemistry model of a helium-air mixture has not been tested on high speed subsonic or 
supersonic on- or off-design jets before. These simulations represent the first steady helium-air mixture jet 
simulations for this flow regime. 

To match experiments, axisymmetric simulations are conducted, using computational grids developed 
and described in more details in a previous study 10 The axisymmetric form of the governing equations are 
solved. Nozzle coordinates are consistent with those used in experiments. Two computational grids are used 
for this study. The first is for the converging nozzle, M d = 1.0 and the second is for the converging- diverging 
nozzle, Md = 1.5. The computational domain extends 75 D from the nozzle exit in the downstream direction, 
50 D from the nozzle exit in the radial direction, and 5 D upstream from the nozzle inlet. Initial iterations 
are performed using a constant CFL number of 0.50 and then explicit Euler time stepping using a constant 
time step is used until the L 2 residual remains constant and a visual inspection of the solution appears 
steady. The L 2 residual is a global measure of the residual error of the governing equations of the flow-field. 
Grid independence studies are performed by using every other grid point in the computational domain in 
both the streamwise and radial directions. No difference in the steady solution between the sequenced (every 
other grid point used) and unsequenced grid is apparent in any of the solutions. The boundary layer at the 
nozzle exit is fully resolved by setting the first grid point in the viscous sublayer region using y + = 1, where 
y + is the coordinate of the viscous wall layer using average conditions. 12 

Boundary conditions for the simulations consist of a downstream outflow in the streamwise direction 
allowing the flow to exit the computational domain. This assumes that the outflow is subsonic and the free- 
held pressure is set as the ambient pressure. Free-stream boundary conditions are specified in the free-held 
where a small uniform how at = 0.001 is assumed for stability and the total ambient pressure is equal 
to the downstream static pressure. A standard axisymmetric boundary condition, which is essentially an 
inviscid wall, lies on the centerline at y/D = 0. Heat transfer is not considered in the study and adiabatic 
no-slip wall boundary conditions are used on the nozzle surface. Finally, the upstream boundary condition 
for the nozzle inlet specihes the total pressure and total temperature of the jet with an initial inlet M = 0.15 
normal to the inlet plane. The inlet M may huctuate and the total values of P and T are held constant. Full 
details of the computational grids and boundary conditions used in this simulation can be found in Miller 
et al. 10 The numerical implementation and mathematical development are shown in Nelson and Power. 20 
These boundary conditions are sufficient for a steady solution to be found when using an ideal gas. Some of 
the simulations use a helium-air mixture, and thus the mass concentration of the helium species is needed 
at the free-stream and inlet boundary conditions. In the free-stream the mass concentration of helium is 
zero and at the nozzle inlet the mass concentration is set to match M a at the nozzle exit. When helium-air 
simulations are performed TTR = 1.0 at the nozzle inlet. 

D. Summary of Jet Operating Conditions 

Experimental measurements were performed in the Pennsylvania State University jet noise facility with pure 
air and helium- air mixture jets. No experimental measurements were performed with heated air jets as the 
facility does not have this capability. CFD simulations were performed using Wind-US for both heated jets 
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and the corresponding helium-air mixtures. The jet conditions experimentally and numerically investigated 
are summarized in Table 1. 

Two nozzles are used for the experimental measurements: a contoured, purely converging (design Mach 
number M ^ = 1.0) nozzle and a converging- diverging nozzle designed with the method of characteristics to 
operate shock- free at a Mach number M& = 1.5. The coordinates of these nozzles are imported directly into 
Gridgen software to create the computational domains for the corresponding Wind-US simulations. 


Table 1. Jet operating conditions for the numerical simulations and experiments. 
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1.26 
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0 

0 

1.4 
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1.5 

1.62 
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1.60 
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44.7 

10.1 

1.49 

467 

363 

1.5 

1.5 

2.22 

637 

1.87 

184 

63.7 

19.5 

1.54 

636 

423 

1.5 

1.5 

3.6 

811 

2.39 

175 

84.1 

42.3 

1.60 

1044 
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III. Results 

This section presents experimental and numerical results from the jet conditions summarized in Table 1. 
The streamwise velocities at two traverse locations, x/D = 4 and x/D = 8, for both the Mj = 0.9 and 
Mj = 1.5 jets at all temperature ratios are shown first, in Figs. 3 and 4. The numerical predictions at 
TTR = 1 are obtained with pure air, while the predictions of the TTR = 2.2 and TTR = 3.6 cases are 
obtained with both heated air and helium air- mixtures. Experimental results are plotted as points, numerical 
results with heated air as solid lines, and numerical results with helium- air mixtures as dashed line. The 
first observation that can be made is that the numerical results from both methods collapse extremely well. 
This illustrates that there is no difference between the predicted velocity profile of a heated air jet and the 
corresponding helium- air jet. This demonstrates that the helium- air mixture provides velocity profiles that 
are similar to the corresponding heated jet condition. The experimentally determined velocity profiles also 
show good agreement with the numerical predictions. However, small discrepancies appear in the inner part 
of the mixing layer, where the experiments generally show lower velocity than the numerical results. CFD 
results of high speed subsonic jets demonstrate lower mixing close to the nozzle exit relative to experimental 
data. In supersonic jets near Mj = 1.5, compressibility effects lower the mixing rate in the experiment, that 
in turn agree better with CFD results. These discrepancies are more visible in the non-dimensionalized plots 
presented later. The centerline velocity obtained from both the CFD and the Pitot measurements increases 
with increasing (simulated) TTR for a jet of given M. These velocities match the theoretical values presented 
in Table 1. In the region of the outer part of the jet, r/D > 0.6, the effects of helium (or heat) addition on 
the velocity profile is much less dramatic: the increase in velocity is very slight, and the curves quickly merge 
into one velocity profile. This is in agreement with measurements performed by Lau 2 with laser velocimetry 
in Mj = 0.5 heated jets of different temperature ratios. 
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Figure 3. Velocity profiles measured in jets of various TTR at x/D — 4. Heated simulations, helium-air 
simulations, and experiments are represented by solid lines, dashed lines, and symbols respectively, a) Mj = 0.90 
b) Mj = 1.50. 




Figure 4. Velocity profiles measured in jets of various TTR at x/D = 8. Heated simulations, helium-air 
simulations, and experiments are represented by solid lines, dashed lines, and symbols respectively, a) Mj = 0.90 
b) Mj = 1.50. 
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It was also found by Lau et al. 15 that the non-dimensionalized velocity profiles for isothermal jets could 
be perfectly collapsed when plotted as a function of the normalized radial distance rj = (r — where 

= Umax/ [dU/dr\ max represents the vorticity thickness and ro.5 is the half velocity point. The collapse 
was shown to follow a Gortler error function profile: 

U 1 /U j = ±[l-erf(ar l )] (12) 

where a is an experimentally determined constant. Similar measurements by Lau 2 showed the same collapse 
for heated jets. Performing a similar non-dimensionalization of both the experimental data and numerical 
predictions of this study, yields a very good collapse of the data, as shown in Figs. 5 and 6. Both the CFD 
and the experimental results collapse very well across all temperature ranges and are therefore self-consistent. 
The helium- air simulation results also match very closely to the heated air results, providing confidence in 
the chemical frozen module of Wind-US. However, this representation highlights the discrepancies between 
the computational results and the experiments. The experimental results fit closely to an error function with 
coefficient a = —2.055 for the Mj = 0.9 jets and a = —2.073 for the Mj = 1.5 jets. On the other hand, the 
numerical results are over-predicting the velocity around r] = ±0.5. This can be attributed to the limitation 
of the Menter SST model to correctly predict the spreading the jet correctly, as previously discussed. 10 

Since the velocity data obtained from experiment can be modeled by an error function, the first and 
second derivative of the velocity profile with respect to radial distance can easily be computed. In the 
previous section, Eq. (7) was solved by making the assumption that the second derivative of the velocity 
profile is negligible compared to the square of its first derivative. The two quantities can now be compared 
by plotting their ratio, as shown in Fig. 7. It is obvious from this plot that for any value r] comprised between 
-0.7 and 0.7 the ratio above mentioned is greater than 10. Since this r] range corresponds to the whole mixing 
layer, this proves that the first term of Eq. (7) can indeed be neglected. The assumption originally made 
that the concentration of helium varies linearly with the streamwise velocity is valid. 

Since Wind-US has the ability to calculate the mass concentration of the helium species of the frozen 
chemistry model spatially, a direct comparison can be conducted with the experiments to further validate the 
assumptions made. The helium mass concentration values are plotted for both experimental and numerical 
data against u m/s in Fig. 8 at radial locations r/D = 4 and r/D = 8. The experimental measurements are 
shown as scatter points and the computation is represented by lines. Clearly, each set of either numerical 
or experimental data shows an almost perfect linear relationship between the velocity and helium mass 
concentration. For the supersonic on- and off-design jets the agreement between the slope of the prediction 
and slope of the experiment is almost perfect. The predictions for TTR = 2.2 and TTR = 3.6 of all three 
Mj, are extremely similar and in the figure can not be differentiated easily. The experimental data agrees 
in a similar fashion except for the Mj = 0.9 cases which show some differences between the x/D = 4 and 
x/D = 8 locations. This is attributed to the fact that the measurements are downstream of the potential 
core, which results in a lower centerline velocity for the experimental data. The predicted potential core 
length does not perfectly match, due to the under-prediction of the spreading rate discussed earlier. As a 
result, the larger centerline velocity creates the slight discrepancy in the slopes observed in Fig. 8. This 
(small) discrepancy is therefore not due to the helium-air frozen chemistry model, which is still believed to 
perform adequately. 
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Figure 5. Non-dimensionalized radial mean velocity profiles measured in fully expanded Mj = 0.9 jets of 
different simulated temperature ratios. 
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Figure 6. Non-dimensionalized radial mean velocity profiles measured in fully expanded Mj = 1.5 jets of 
different simulated temperature ratios. 
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Figure 8. Comparisons of the simulated and experimental helium mass concentration at x/D = 4 and x/D = 8 
as a function of u m/s. Heated simulations, helium-air simulations, and experiments are represented by solid 
lines, dashed lines, and symbols respectively. 
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The helium concentration can also be plotted as a function of the radial position, as shown in Fig. 9 for 
the experimental conditions of Table 1 and the corresponding helium- air mixture simulations. The values of 
the helium species are shown at the two downstream measurement locations of x/D = 4 and x/D = 8 . Good 
overall agreement between the predictions and experiment can be seen on the centerline and at approximately 
r/D = 0.5 for all the jet conditions examined. Unfortunately, as mentioned before, the Menter SST model 
does not correctly predict the jet spreading rate, resulting in some discrepancies very similar to the ones 
observed in Figs. 3 and 4. These discrepancies are not due to the additional use of the frozen chemistry 
model for helium-air mixtures. Unfortunately, the steepened meanflow profiles of the predictions relative to 
the experimental results has an effect on the helium mass concentration. 



M=0.9 TTR=2.2 



,, M=0.9TTH=2.2 


Figure 9. Comparisons of the simulated and experimental helium mass concentration as a function of r/D. 
a) x/D — 4 b) x/D = 8. Helium-air simulations and experiments are represented by dashed lines and symbols 
respectively. 


Cross-stream temperature profiles of the Mj = 0.9 jets can be seen in Fig. 10. The solid lines correspond 
to the static temperature, as obtained from heated air numerical simulations. The experimental points 
correspond to the static temperature as simulated by the local helium concentration. Similar to the profiles 
of Ui, the spatial variables have been normalized by the nozzle diameter D. Dimensional units of Kelvin have 
not been normalized in order to illustrate the differences in simulated temperatures of the various Mj =0.9 
jets. Only the cross-stream profiles at x/D = 4 are shown, since the simulations become less reliable further 
downstream. The agreement is excellent at all TTR. Since the simulations were performed for a heated air 
jet, this constitutes a very good validation of the processing approach used for the experimental data as well 
as the ability to properly simulate heat with the addition of helium. 

If the temperature profiles are normalized by Tj and plotted against 77 then self-similarity can be shown 
for the temperature like the velocity. This is shown in Fig. 11 for the Mj = 0.9 jet at TTR = 2.2 and 
TTR = 3.6. At both TTR the collapse of the numerical predictions and experimental data is very good. 
Note that 77 is based on S w and not the maximum temperature gradient in the cross-stream direction. This 
fact causes the temperature profiles to not have values of 0.5 at 77 = 0.0 like the velocity plots. 
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Figure 10. Temperature profiles measured in the Mj = 0.9 jets at x/D = 4. Heated simulations and experiments 
are represented by solid lines and symbols respectively. 
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Figure 11. Non-dimensionalized radial mean temperature profiles measured in the fully expanded Mj = 0.9 jet 
of different simulated temperature ratios. The lines are simulated heated air and the symbols are experimental 
data. 
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Since the local speed of sound is calculated at each radial location of the experimental data, local Mach 
number (Mi) profiles can be compared with corresponding numerical simulations, as shown in Fig. 12. 
Figure 12 shows the radial distribution of Mi at x/D = 4 in part a) and x/D = 8 in part b). Both the 
simulated heated air and helium-air mixture results are also shown. As seen in the cross-stream plots of 
E/i, the comparison between the numerical and experimental results shows slight over prediction of Mi on 
the inner and outer parts of the mixing layer, similar to the result for the local velocity U\. Slight over- 
prediction of Mi on the inner and outer parts of the mixing layer are apparent, as it is for the local velocity 
U\. Similarly, the jet centerline Mach number is obviously over-predicted at x/D = 8. This is all consistent 
with the fact that the Menter SST model under-estimates the growth of the mixing layer. The excellent 
agreement between helium-air and heated air simulations again brings confidence in the accuracy of the 
frozen chemistry model of Wind-US and the ability to properly simulate the flow of a heated jet with the 
addition of helium. 



0.8 
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0,2 
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Figure 12. Mach number profiles of the numerical and experimental data for the M = 0.90 jet operating at 
various temperature ratios, a) Radial location x/D = 4. b) Radial location x/D = 8. Heated simulations, 
helium-air simulations, and experiments are represented by solid lines, dashed lines, and symbols respectively. 

Schlieren images of the Mj = 1.5 jet operating on- and off-design are obtained using helium-air mixtures 
to simulate the various TTR. Corresponding images from numerical prediction can be constructed by post 
processing the density flow-field. Constructing these numerical Schlieren requires the derivative of p in the 
cross-stream direction to be calculated. Since the simulations are axisymmetric the result is mirrored across 
the x-axis. In a cartesian x — y coordinate system, because the values of dp/dy are essentially opposite 
across the x-axis for an axisymmetric jet, the values of dp/dy are multiplied by negative one in the lower 
half plane. This creates the same light and dark gradients as observed from the Schlieren measurements. 
The contours of the numerical results are continuous and range from constant values of -6 to 6, though the 
values of the density gradients are much larger. These values are selected to better match the colors of the 
experimental plots. Numerical Schlieren obtained show very good overall agreement with the experimental 
axisymmetric jets. If the jet is not axisymmetric then three-dimensional flow features may not be apparent 
in the numerical Schlieren due to the fact that only a computational plane of data may be utilized to create 
the image (while the experimental Schlieren are an integrated representation of the whole jet). 

A comparison of the supersonic M^ = 1.0, Mj = 1.5, TTR = 2.2 jet is shown in Fig. 13(a). The top 
image in a) is produced experimentally with a Z-type Schlieren setup described in Veltin and McLaughlin. 21 
The images presented result from a large number of spark Schlieren images averaged in order to obtain a 
better visualization of the shock cell structure and spreading angle. As a side effect, small scale structures 
are less apparent. Averaging is not required to produce numerical Schlieren. The middle and bottom image 
of part a) is representative of the numerical Schlieren. The middle numerical Schlieren uses heated air 
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and the bottom image shows the helium-air jet for a simulated TTR = 2.2. Qualitative comparisons may 
be made by examining the images. In general the on-design case shows good overall agreement with the 
experiment. The shock cell structure may be faintly seen in the experimental Schlieren but are much more 
apparent in the numerical simulation. Wind-US inherently assumes that the boundary layer is turbulent 
and never laminar and no-tripping of the boundary layer occurred in the nozzle interior. This is obviously 
not representative of the physics of these small jets with Reynolds number around 400,000 that probably 
possess laminar boundary layers within the entire nozzle length. Since turbulent boundary layers have a 
larger displacement thickness than laminar ones, the simulation produces slight shocks when operating at 
Mj = 1.5. The convergent-divergent nozzle used in the experiment is designed to operate with a fully 
expanded jet and a laminar boundary layer at M ^ = 1.5. The shock cell length also exhibits a slight change 
when operating with a helium-air mixture. Both jets are very close to Mj but are slightly different due 
to very small variations in the flow properties (such as a slightly decreased viscosity due to a different gas 
species in the nozzle). The off-design parameter of the shock-cell structure is extremely small. 

A second Schlieren comparison is made in Fig. 13(b) and illustrates the heated off-design case when 
Md = 1.0, Mj = 1.47, and TTR = 2.2. The arrangement of the images is the same as part a). This case 
uses the convergent conical nozzle in both the Schlieren experiment and simulations. The same method to 
produce the numerical Schlieren of Fig. 13(a) is also used here with the same continuous contours and range. 
Excellent qualitative agreement may be observed with regard to the shock cell structure position between 
the predictions and experimental results. 



(a) M d = 1.50, Mj = 1.50, TTR = 2.20 jet. (b) M d = 1.00, Mj = 1.50, TTR = 2.20 jet. 


Figure 13. Comparisons of Z-type Schlieren with numerically generated Schlieren. The top images are pro- 
duced by the Z-type experimental Schlieren and the middle and bottom are numerically generated results 
produced by heated air or helium-air mixtures. 
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Mach number contour plots are created from the numerical results of the Mj =0.9 and Mj = 1.5 heated 
jets at TTR = 2.2 and TTR = 3.6. These are illustrated in Figs. 14 through 17. In each figure the top 
half plane is representative of a simulated heated air result and the bottom half plane is the helium-air 
result. For the Mj =0.9, TTR = 2.2, and TTR = 3.6 jet the qualitative agreement between the two jets 
is excellent. The core lengths and shear layer growth of both the heated air and helium- air mixture are 
almost the same. The same good agreement of core length and shear layer growth may also be seen in the 
supersonic cases of Figs. 16 and 17, however, there is some discrepancy between the shock cell structures 
similar to the observation made from the Schlieren images. These once again, can be attributed to thicker 
boundary layers and changes in 7 due to helium in contrast to air. 



Figure 14. Contour plot of M for the M ^ = 1.0, Mj = 0.90, TTR = 2.20 jet. The top half plane, y > 0 are contours 
of the heated simulated air jet while the bottom half plane, y < 0 is representative of contours of the helium-air 
jet. 



Figure 15. Contour plot of M for the M ^ = 1.0, Mj = 0.90, TTR = 3.60 jet. The top half plane, y > 0 are contours 
of the heated simulated air jet while the bottom half plane, y < 0 is representative of contours of the helium-air 
jet. 
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Figure 16. Contour plot of M for the M d = 1.50, Mj = 1.50, TTR = 2.20 jet. The top half plane, y > 0 are 
contours of the heated simulated air jet while the bottom half plane, y < 0 is representative of contours of the 
helium-air jet. 



Figure 17. Contour plot of M for the M ^ = 1.50, Mj = 1.50, TTR = 3.60 jet. The top half plane, y > 0 are 
contours of the heated simulated air jet while the bottom half plane, y < 0 is representative of contours of the 
helium-air jet. 
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Centerline plots are produced for pi, U\, and M\ from simulations of Mj = 0.9 and Mj = 1.5 heated and 
heat-simulated jets at TTR = 2.2 and TTR = 3.6 in Figs. 18 and 19 respectively. In both figures the rows 
of sub-figures are representative of pi, Ui, and Mi while the columns correspond to the different operating 
conditions. Since the helium- air concentrations are chosen to match the exit acoustic Mach number M a 
for both TTR values, the values of M\ at the exit of the nozzle (x/D = 0) matches perfectly between the 
mixture and the heated air simulations. For example, in part g) of Fig. 18 M\ is almost identical along the 
centerline axis. The same excellent agreement can be seen in other cases also. Only for the supersonic cases 
does the core extend slightly farther with the helium-air jets compared to the heated air jets. As expected, 
comparisons of the density for both temperature ratios and three corresponding jet conditions do not exhibit 
the same level of agreement. This is due to the fact previously stated that it is not possible to match both 
the density and the Mach number. While these simulations are performed by matching M a between heated 
air and helium-air jets, the calculation could easily have been performed by matching the exit densities. 
Nevertheless, for each TTR and NPR the discrepancy in the predicted values of pi remains very small. 

The values of the streamwise velocity, C/i, on the centerline also agree very favorably, as seen in Figs. 18 
and 19 for both TTR and Mj values. These plots in particular illustrate the very slight changes in core 
length of the off-design calculation. Still, the core length remains almost exactly the same between the 
heated air and helium-air mixture for the subsonic calculations and supersonic on-design calculations. 

The centerline values of the supersonic on-design jet of heated air and helium- air demonstrate that the 
core lengths are exactly the same. Also, it can be seen in part h) of both figures that the values of Mi at the 
jet exit are almost exactly the same value. Even with both of these observations, for the supersonic on-design 
jet, the helium-air jet is slightly over-expanded while the heated air simulated jet is under-expanded. The 
rise or fall of velocity due to shocks or Prandtl-Meyer expansion waves can be seen in the centerline plots. 
These values of Mj at the jet exit are extremely close to the design Mach number M ^ = 1.5. Their slight 
difference produces shock cells of different lengths, as already observed from the Schlieren images and the 
Mach number contour plots. When examining the discontinuities of Mi in part h) of both Figs. 18 and 19, 
it appears that there is an initial very strong shock. However, this is an artifact of using the axisymmetric 
boundary condition where dPi/dr = 0. This condition is theoretically true at a point; however, it is enforced 
in CFD by using grid points that are on the centerline axis and the corresponding grid points of the spatial 
discretization stencil in the cross-stream direction. Therefore, there is a numerical effect of broadening the 
region where dPi/dr = 0. This effect can be varied by changing the grid spacing in the cross-stream direction 
near the centerline, effectively causing a barrel shock that has the diameter on the order of twice the distance 
from the centerline to the first grid point in the cross-stream direction. To minimize this effect, the grid 
point clustering in the cross-stream direction on the centerline axis uses the same distribution as the nozzle 
wall at x/D = 0. 
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Figure 18. Various centerline comparisons of TTR = 2.2 jets. air helium-air. 
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Figure 19. Various centerline comparisons of TTR = 3.6 jets. air helium-air. 
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IV. Conclusion 


Experimental and numerical assessment of the flow properties of high speed helium-air mixture and heated 
air jets were presented. High subsonic and supersonic shock-free and shock- containing jets were considered 
and demonstrated a very good agreement between the numerical simulations and experiments. The flow 
experimentally assessed in helium-air jets has properties very similar to heated air jets measurements present 
in the literature. The turbulent model used for the simulation was observed to underpredict the growth rate 
of the jet in particular for the lower jet Mach number cases, leading to a mis-prediction of the potential core 
length. However, the very good agreement between simulated heated air and helium-air mixture jets provide 
a validation of the frozen chemistry model used, as well as further evidence that the helium addition does 
indeed provide mean flow properties similar to heated air. 
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